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ABSTRACT 


Several 3-D finite element analyses of (O/0/-0) g graphite epoxy laminates, 

where 0-15, 20, 25, 30, and 45 degrees, subjected to axial tension load, were 
performed. The interlaminar stresses in the 0/-0 interface were calculated 
with and without a matrix crack in the central -0 plies. The interlaminar 
normal stress changes from a small compressive stress when no matrix crack is 
present to a high tensile stress at the intersection of the matrix crack and 
the free edge. The analysis of local delamination from the -0 matrix crack 
indicates a high strain energy release rate and a localized mode I component 
near the free edge, within one ply distance from the matrix crack. In order 
to examine the stress state causing the matrix cracking, the maximum 
principal normal stress in a plane perpendicular to the fiber direction in 
the -0 ply was calculated in an uncracked laminate. The corresponding shear 
stress parallel to the fiber was also calculated. The principal normal stress 
at the laminate edge increases through the ply thickness and reached a very 
high tensile value at the 0/-0 interface indicating that the crack in the -0 
ply may initiate at the 0/-0 interface. Crack profiles on the laminate edge 
in the -0 ply were constructed from the principal stress directions. The 
cracks were found to be more curved for layups with smaller 0 angles, which 
is consistent with experimental observations in the literature. 


1 . INTRODUCTION 


In laminated composite materials, matrix cracking and local delaminations 
from matrix cracks are important damage mechanisms which contribute to 
stiffness loss and eventual laminate failure [1,2]. Therefore, it is 
Important to understand and predict the Initiation and growth of this damage. 

In the static tension test of a [0^/30 4 /-30^] g graphite epoxy laminate, Soni 

and Kim [3] observed a curved crack in the -30 degree group of plies. The 
tests conducted on [0 /±15] g laminates by Lagace and Brewer [4] showed that 

the delamination area in the 15/-15 interface extended between the transverse 
ply crack and the free edge. O'Brien and Hooper [5] found that in (Oj/^/" 

0 y) graphite epoxy laminates (0=20, 25, 30 degrees) subjected to tension 

fatigue loading, matrix cracking was detected as the first event followed by 
local delamination. The matrix cracks formed near the stress -free edge and 
the local delaminations formed in the 0/-0 interface. These local 
delaminations were bounded between the free edge and the matrix crack (figure 
1) . In an experimental and analytical study, Fish and O'Brien [6] concluded 
that in (+15/90^/- 15) g glass epoxy laminates the interlaminar tensile 

stresses due to cracking in the +15 degree ply plays a significant role in 
the onset of local delaminations. 

The purpose of this paper is to study the stress state causing matrix 
cracking and delamination in the (0 /0/-8) graphite epoxy laminates subjected 

to axial tension. Three dimensional finite element analyses were performed on 
(O/0/-O) laminates for 0=15, 20, 25, 30, and 45 degrees. The effect of the 
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- 8 degree ply crack on the interlaminar stress in the 8/ -8 interface was 
evaluated. A strain energy release rate analysis of delamination from the - 6 
ply crack was performed. In order to better understand the matrix cracking, 
the in-plane stresses in the - 6 ply of the uncracked laminate were 
calculated. The principal stresses in the -6 ply were examined to predict the 
probable crack profile. The crack profiles predicted using the present 
analysis were compared with those from the experiments in the literature [5], 

2. LAMINATE CONFIGURATION AND LOADING 

A (0 / 8 /- 0 ) graphite epoxy laminate subjected to axial tension load (figure 

2) was analyzed using the finite element method. The laminate length was 
assumed to be 75h and the width was lOh, where h is the ply thickness. A unit 
ply thickness was assumed in the model. The Cartesian coordinate system has 
its origin at point A and the axes are shown In figure 2. The matrix crack 
was assumed to occur in the two - 8 degree central plies, as evidenced in the 
experiments [5]. The direction of the matrix crack across the laminate width 
was assumed to be parallel to the fiber direction. For simplicity, the 
matrix crack was assumed to be straight through the thickness of the - 8 ply 
and oriented normal to the XY plane. The laminate was assumed to be 
constrained at x=0 plane and an arbitrary uniform axial displacement of 
0.5625h was applied on the x=75h face which corresponds to a strain of 
0.0075. 

The finite element discretization of the laminate used in the analysis is 
shown in figure 3. Only the top half of the laminate (above z=0) was modelled 
and appropriate symmetry boundary conditions were applied on the z-=0 plane. 
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In order to analyze the effect of the stress free edge on matrix cracking and 
on local delamination, a finely graded mesh was used at one of the free edges 
(y=0) . The smallest element size used in each direction was .Olh and the mesh 
refinement is assumed to be adequate for stress analysis accuracy. The 
behavior of the matrix crack in the (0 / 8 /- 9 )^ laminate and the resulting 

local delaminations depend on the angle 8 of the off-axis plies [5] . 
Therefore, the stress analysis was carried out for each of the laminates, 
where the angle 8 was assumed to be 15, 20, 25, 30, and 45 degrees, A 
separate model was constructed for each case by assuming a matrix crack plane 
at the corresponding angle - 8 (plane BCDE in figure 2) and a finely graded 
mesh was used on either side of the plane BCDE. The 8 and - 8 degree plies had 
six elements each through the thickness, and the mesh was finely graded near 
the 8/ -8 interface. A total of 4312, 20-node isoparametric tetrahedral 
elements having 20477 nodes were used in each model. The matrix crack and 
delamination were simulated by releasing coincident multi-point constraints 
provided in the model. The elastic constants for the graphite fiber and epoxy 
resin matrix material used in the analysis are shown in Table 1 [5]. 

3. ANALYSIS 

Several 3-D finite element analyses were performed in order to study various 
aspects of matrix cracking and local delaminations initiating from the matrix 
cracks. The various analyses are described below. 

3.1 INTERLAMINAR STRESSES : The (0/15/-15) laminate was tested in 

reference 4 and the delamination was attributed to the interlaminar shear. 
Therefore the 8 ~ 15 degree case was analyzed here. 
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The experimental evidence [5] indicated that the local delamination on the 
0/-0 interface was associated with matrix cracking in the -0 degree ply. The 
crack in the -15 ply was assumed to be straight in the thickness direction, 
oriented normal to the XY plane, and parallel to the fiber direction. The 
interlaminar normal and shear stresses were obtained in the 15/-15 interface 
in the laminate before and after introducing a matrix crack in the -15 ply, 
to examine any change in these stresses due to the introduction of the matrix 
crack, 

3.2 STRAIN ENERGY RELEASE RATE: Local delamination in the 45/-45 
interface was modelled as initiating from a matrix crack in the -45 ply. The 
strain energy release rate (G) was calculated for small increments of 
delamination growth. The delamination front was assumed to be uniform and 
parallel to the matrix crack for convenience of modelling. The elements were 
in the shape of a parallelopiped, one side being parallel to the free edge. 
The virtual crack closure technique (VCCT) was used to compute the strain 
energy release rate [7]. The equations used in the G calculation are also 
given in Reference 8 . 

3.3 IN-PLANE STRESSES: The normal tension stress on the BCDE plane In an 
uncracked laminate was examined for 0=15, 30, and 45 degrees. If the tension 
stress extends only partially along the potential crack plane, then the 
matrix crack may extend only partially across the laminate width. 

3.4 PRINCIPAL TENSION AND ACCOMPANYING SHEAR: The curved shape of the 
matrix crack on the stress free edge (y-0) in the -0 ply, observed in the 
experiments [5], may be due to the principal tension stress, rather than just 
the in-plane tension and In-plane shear stresses. The principal tension 
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stress may be acting at an angle other than 90 degrees to the vertical plane 
BCDE (figure 2). The principal tension stress acting perpendicular to the 
fiber was computed as explained below. 

The 3-D finite element analysis of the intact laminate was performed using 
the appropriate mesh for the ply angle (0) under consideration. The stresses 
at the nodes on the potential matrix crack surface BCDE (figure 2) were then 
transformed to the coordinate system 123 (figure 4) on plane FGHJ which is 
normal to plane BCDE. In the system 123, axis 1 was normal to the vertical 
plane BCDE, axis 2 was parallel to BCDE and the axis 3 was vertical (parallel 
to z direction). The principal normal stress due to the three stress 
components (a^, the plane were computed at each node. The 

corresponding direction of the principal stress, i.e. , the angle a measured 
counter-clockwise from axis 1, was also determined at each node. The six 
stress components in the 123 coordinate system were then transformed through 
the angle a in the 1'3' plane so that axis 1' coincides with the principal 
direction. 

The new coordinate system 1'2'3' is also shown in figure 4. The purpose of 
this second transformation was twofold. The normal stress on t * le 

transformed plane in the 1' direction should check with the principal stress 
calculated earlier. In addition, the shear stress acting parallel to the 
fiber direction on this plane was obtained. Thus, the principal 

tension stress normal to the fiber direction, the angle a indicating the 
direction of the principal stress in the 1'2' plane, and the shear stress 
acting parallel to the fiber direction were obtained for each node on plane 
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BCDE. This analysis was repeated for the five types of laminates, with ply 
angles, 0, of 15, 20, 25, 30, and 45 degrees. 

4. RESULTS AND DISCUSSION 

4.1 Interlaminar Stresses: The interlaminar normal stress, o . in the 

z 

+15/- 15 interface, normalized by the remote applied axial stress, was 

found to have a small compressive value near the free edge of the (0/15/- 15) 

uncracked laminate [5]. The 3D analysis of the laminate with a matrix crack 
in the -15 degree ply, however, indicated that there was a large tensile 
interlaminar normal stress at the intersection of the matrix crack and the 
stress free edge (fig. 5). This high interlaminar tension will contribute 
significantly to the onset of delamination in the 15/-15 interface (figure 
1 ). 

The interlaminar shear stress in the +15/-15 Interface of the uncracked 

laminate increases in magnitude towards the stress free edge but is uniform 
along any line parallel to the free edge. After modelling a crack in the -15 
degree ply, the interlaminar shear stress increased in the area bounded by 
the matrix crack and the free edge (figure 6). The increase in interlaminar 
shear stress was less than one order of magnitude. 

4.2 STRAIN ENERGY RELEASE RATE: A uniform delamination front growing 

parallel to a -45 degree ply crack in a (0/45/-45) s laminate was modelled. In 

this analysis, the direction of delamination growth is not normal to the 
delamination front and the delamination growth is modelled between dissimilar 
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materials. Therefore, this modelling was expected to give only a qualitative 
indication of the effect of the free edge on the strain energy release rate. 

2 

The variation of mode I component of G and the total G normalized by e tE^ 

with distance s along the delamination front, for a/h-0.02 are shown in 
figure 7, The symbols e, t, and denote the remote axial strain, the 

laminate thickness and Young's modulus in the fiber direction, respectively. 
Both the total G and G^ distributions show a marked increase within one ply 

thickness of the free edge. The Mode I component appears to be absent in the 
interior along the laminate width. Figure 8 shows the plots of normalized G 
versus a/h near the free edge of the laminate and in the interior of the 
laminate. G decreases as the delamination grows away from the matrix c rack. 
Figure 8 also indicates that G is higher near the free edge than in the 
interior. 

4.3 IN - PLANE STRESS: The plane BCDE parallel to the fiber direction shown 

in figure 2 is a potential crack plane in the -0 degree ply. The stress 
component normal to the plane BCDE was obtained for each point (node) on this 
plane. The points on plane BCDE where the normal stress changes from tension 
to compression are marked and joined to indicate the area under in-plane 
tension. The projection of plane BCDE and its area under tension, on yz plane 
is shown in figure 9, for laminates with 0*15, and 30 degrees. In order to 
obtain the area of plane BCDE under tension along the crack which is parallel 
to the fiber direction, the y coordinate must be increased by a factor of 
l/sin(0). Figure 9 shows that for 0*15, and 30 degrees the in-plane tension 
exists up to a distance of 1 to 1.5 plies in the y direction on yz projection 
or three to four plies in the -0 degree ply. This result is consistent with 
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the finding from quasi-3D analysis in ref. 5. Furthermore, figure 9 shows that 
this in-plane tension extends through the thickness to the 9/ -9 interface 
only near the free edge. Therefore, the delamination may initiate between the 
9 and -9 plies where - 9 ply meets the laminate edge. The in-plane normal 
stress for the (0/45/-45) s laminate, however, was tensile over the entire 

crack face. Hence, the matrix crack in the -45 ply and the delamination in 
the 45/-45 interface from -45 ply crack may extend further towards the 
interior of the laminate than for 0-15 - 30 degree cases. 

4.4 MATRIX CRACK PROFILE: The angles (a) with respect to horizontal axis 
1 (figure 4), at which the principal stress acts near the laminate edge were 
calculated for fiber angles 15,20,25,30, and 45 degrees. The values of a 
increase from the center of the laminate towards the 9/- 9 interface. Thus the 
plane on which the principal normal stress acts tends to move away from the 
vertical axis as z changes from 0 to 1 . By using the values of a for various 
z coordinate values, it is possible to construct an approximate shape of the 
matrix crack caused by the principal normal tension stress The crack 

shapes constructed using this approach are shown in figure 10(a). The cracks 
appear to be curved towards the 9/ -9 interface. Figure 10(a) indicates that 
the crack shape for 0=45 degrees tends to be less curved than that for the 
smaller 9 angles considered. 

The crack profiles in figure 10(a) were constructed on FGHJ planes normal to 
the fiber direction (fig. 2). In order to obtain the shape of the crack on the 
free edge, the c rack must be projected using the corresponding fiber angle. 
This is achieved by multiplying the abscissa values in figure 10(a) by a 
factor of (l/sin0). The crack profiles on the free edge thus obtained are 
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shown in figure 10(b). Figure 10(b) clearly shows that as the angle 0 
decreases, the crack profile tends to be more curved. These analytical 
results exhibit the same trend as seen in the micrographs of the specimen 
edges from experiments on similar laminates [5], and are sketched in figure 
11 for clarity. 

During the course of the present analysis, the authors were made aware of 
the paper on the matrix cracking in cross -ply laminates by Groves et. al. 
[9]. Reference 9 also used the idea of principal stresses to address the 
occurrence of curved cracks in 90 degree plies and observed a similar 
correlation between the crack profile and direction of principal stress. 

4.5 PRINCIPAL TENSION STRESS AND SHEAR STRESS: The values of the 
principal tension stress acting in a plane normal to the fiber direction at 
the free edge were calculated for fiber angles 15,20,25,30, and 45 degrees. 
The values increase from the laminate midplane (z=0) towards the 0/-0 

interface (z=l) for each of the laminates analyzed. The shear stress along 
the fiber direction, r increased from the laminate midplane towards the 

0/-0 Interface, through the ply thickness for all values of 0 as well. 

The principal tensile stress and the shear stress calculated above were 
plotted in figure 12 normal to the expected crack profiles shown previously 
in figure 10(a). As seen in figure 12, in the case of 0=15 degrees, the value 
of is smaller than the shear stress whereas the two are almost equal 

for 0=45 degrees. For the laminates considered, the va ^ ue near t ^ xe 

interface appears to be highest (0.24 GPa) for 0=25 degrees. 
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The results discussed above indicate that the principal tension stress normal 
to the fiber and the shear stress parallel to the fiber near the laminate 
edge are highest at the interface. Therefore, the matrix crack is most likely 
to form in the -8 ply near the Q/-8 interface. 

5. CONCLUDING REMARKS 

A (O/0/-0) g graphite epoxy laminate subjected to axial tension was analyzed 

using three dimensional finite element analyses. The first analysis was 
performed to examine the effect of matrix cracking in the -15 degree ply on 
the interlaminar stresses between the 15/-15 interface of a (0/15/-15) s 

laminate, particularly near the stress free laminate edge. Next, a strain 
energy release rate analysis was performed for delamination growth from a 
matrix crack. Then, the stress state in the uncracked laminate was analyzed 
to understand its effect on the crack profile in the -B ply. Based on these 
analyses, the following conclusions were reached. 

The interlaminar normal stress, a , changed in the 15/-15 interface of a 

(0/15/- 15) laminate from a small compression stress when no matrix crack was 

present to a large tension stress near the free edge after the formation of 
cracking in the -15 degree ply. This indicates that the delamination from the 
matrix crack may initiate at the point where the matrix crack intersects the 
laminate edge. 

For the 8 = 45 degree case, the strain energy release rate for a local 
delamination growing uniformly away from the matrix crack in the -45 degree 
ply is higher near the laminate edge than in the interior of the laminate. 
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The value of total G decreases as the delamination grows away from the matrix 
crack. The mode I component of G exists only near the free edge within one 
ply distance from the matrix crack. 

The in-plane stresses in the -0 ply for the case of 0-15, and 30 degrees, in 
the uncracked laminate are tensile only up to a distance of 3 to 4 plies from 
the laminate edge along a potential crack face. The area of tension stress 
extends through the -8 ply thickness, to the 0/-0 interface, only near the 
free edge. In contrast, the in-plane tension stress In the -0 ply for the 
0=45 degree case extends across the ply thickness to the 0/-0 interface and 
throughout the laminate width. This suggests that the cracking in the -45 
degree ply may extend further across the laminate width. 

The principal tension stress was calculated In a plane normal to the fiber 
direction at the laminate edge. This principal stress Increased along the 
thickness of the -0 ply, from the center of the laminate towards the 0/~0 
interface, where it increased sharply. The shear stress acting along the 
fiber direction shows similar variations. Therefore, matrix cracking Is 
likely to initiate at the free edge, near the 0/-0 interface. 

The principal tension direction is inclined at an angle to the plane of the 
laminate. This angle increases from the center of the laminate along the -0 
ply thickness, towards the 0/-0 interface. The projection of the matrix crack 
on the free edge indicates a larger curvature for smaller ply angles. This 
observation is consistent with the experimental findings in the literature. 

To summarize, the occurrence of a high Interlaminar tensile normal stress in 
the 0/-0 interface due to matrix cracking in the -0 ply and the high value of 
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G which is localized in the vicinity of the free edge and matrix crack, will 
result in delamination onset in the 8/ -8 interface. The matrix cracking in 
the - 8 ply appears to initiate at the 0/-8 interface due to high tensile 
principal stress, resulting in curved crack profiles. These crack geometries 
must be modelled to accurately determine interlaminar stresses and G's 
associated with local delaminations. These solutions may then be used in 
progressive failure methodology for laminated composites. 
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TABLE 1 


MATERIAL PROPERTIES 
AS4/3501-6 GRAPHITE EPOXY 


“'ll 


E 22 ’ 

E 33 

G 12 ’ 

G 13 

G 23 


l/ 12’ 

"13 

v 23 



- 134.45 GPa, (19.5 x 10 6 psi) 

- 11.03 GPa, (1.60 x 10 6 psi) 

- 5.84 GPa, (0.847 x 10 6 psi) 

- 2.98 GPa, (0.432 x 10 6 psi) 

- 0.301 

- 0.490 
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Fig. 2- Configuration and loading of a (0 /6/-6) s graphite epoxy laminate 
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H E . 5- Normalized interlaminar normal stresa in 15/-15 interface due to 
matrix crack in -15 ply of (0/15/-15) g laminate 
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Fig. 6- Normalized interlaminar shear stress in 15/-15 interface due to 
matrix crack in -15 ply of (0/15/-15) s laminate 
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Fig. 7- Variation of G across delaraination front in 4S/-45 interface of 
(0/45/- 45 ) s graphite epoxy laminate with a -45° ply crack 
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Fig. 8- Change in G due to delamination growth in (0/45/-45) s graphite 







(a) Crack profile normal (b) Crack profile on 
to crack plane free edge 


Fig. 10- Predicted crack profile in -d plies of (0/9/-6) s laminate 
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laminates following cyclic loading (Ref. 



0.0075 
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stress acting along the fiber direction in -6 ply of (0 /B/-6) laminates 




1. Report No. 

NASA TM- 104062 
A V.QfYYM TR-Q1 -R-002 

2. Government Accession No 

4. Title and Subtitle 

Analysis of Matrix Cracking and Local Delamination in 
(0/0/— ©)s Graphite Epoxy Laminates Under Tension Loading 

7. Author(s) 

S. A. Salpekar* and T. K. O'Brien** 

9. Performing Organization Name and Addre 

NASA Langley Research Center, 
U.S. Army Aviation Research Tec 
Aerostructures Directorate 
Hampton. VA 23665-5225 

ss 

Hampton, VA 23665-5225 
:hnology Activity (AVSCOM) 

idministration 

nmand 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space A 
Washington, DC 20546 
U.S. Army Aviation Systems Coi 
St. Louis, MQ 63166 


NASA 

•AtU'xul Armpit*. ■ . »rxJ 

',1* ilk * 


Report Documentation Page 


3. Recipient's Catalog No. 


5. Report Date 

March 1991 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

505-63-01-05 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

*S. A. Salpekar, Analytical Services and Materials, Inc., Hampton, VA 23666 
**T. K. O'Brien, U.S. Army Aerostructures Directorate, USAARTA-AVSCOM, Langley Research 
Center, Hampton, VA 23665 


16. Abstract 

Several 3-D finite element analyses of (O/0/-9) s graphite epoxy laminates, where 0 = 15, 20, 25, 30, and 
45 degrees, subjected to axial tension load, were performed. The interlaminar stresses in the 0/-0 

interface were calculated with and without a matrix crack in the central -0 plies. The interlaminar normal 
stress changes from a small compressive stress when no matrix crack is present to a high tensile stress at 

the intersection of the matrix crack and the free edge. The analysis of local delamination from the -0 
matrix crack indicates a high strain energy release rate and a localized mode I component near the free edge 
within one ply distance from the matrix crack. In order to examine the stress state causing the matrix 

cracking, the maximum principal normal stress in a plane perpendicular to the fiber direction in the -9 ply 
was calculated in an uncracked laminate. The corresponding shear stress parallel to the fiber was also 
calculated. The principal normal stress at the laminate edge increases through the ply thickness and reachecj 

a very high tensile value at the 0/-0 interface indicating that the crack in the -0 ply may initiate at the 
0/-0 interface. Crack profiles on the laminate edge in the -0 ply were constructed from the principal stres<| 

directions The cracks were found to be more curved for layups with smaller 0 angles, which is 



Angle ply 
Matrix crack 
Delamination 

Three-dimensional finite element analyses 


Unclassified - Unlimited 
Subject Category - 24 


rai 

Oe 


onenen; 


lecurity Cl' 

Unclassified 


rate 

►if.loT this report) 


20 Security Classif. (of this page) 

Unclassified 


["21 No. of pages 
28 


| 22. Price 

AO 3 


NASA FORM 1«26 OCT 86 




